EPOXI at Comet Hartley 2 
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Abstract 

Understanding how comets work, i.e., what drives their activity, is crucial to using 
comets to study the early solar system. EPOXI flew past comet 103 P/Hartley 2, one with 
an unusually small but very active nucleus, taking both images and spectra. Unlike large, 
relatively inactive nuclei, this nucleus is outgassmg primarily due to CO;, which drags 
chunks of ice out of the nucleus. It also shows significant differences in the relative 
abundance of volatiles from various parts of the nucleus. 
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introduction 

Comets are the fundamental building blocks of the giant planets and may be an 
important source of water and organics on Earth. On 4 July 2005, the Deep Impact 
mission carried out an impact experiment on comet 9P/Tempel 1 (7,2) to study 
differences between the comet’s surface and the interior. Although the impactor 
spacecraft was destroyed, the flyby spacecraft and its instruments remained healthy in its 
3-year, heliocentric orbit after completion of the mission. The Deep Impact flyby 
spacecraft was retargeted to comet 103P/Hartley 2 as part of an extended mission named 
EPOX1. 

The flyby spacecraft carries the High Resolution Instrument (HR1), which combines a 
visible-wavelength camera having pixel size 2 prad and a set of filters with a near-IR 
(1 .05-4.85 pm) spectrometer having an entrance slit 10 prad by 256 mrad with 512 
spatial pixels along the slit. Spectral maps are created by scanning the slit across the 
comet while taking a sequence of spectra. The Medium Resolution Instrument (MRI) has 
a pixel size 10 prad and a different, but overlapping set of visible- wavelength filters (3,4). 


Encounter with Hartley 2 

Closest approach to Hartley 2 was 694 km at 13:59:47.3 1 UTC on 4 November 2010. 
one week after perihelion passage and at 1.064 AU from the sun. Flyby speed was 12.3 
km/s and the spacecraft flew under the comet with a somewhat northward trajectory in a 
solar system reference frame. Because instruments are body-mounted on the spacecraft, 
the spacecraft rotated to keep the instruments pointed at the comet. Observations of the 
comet were carried out for 2 months on approach (5 Sep - 4 Nov) and for 3 weeks on 
departure (4 - 26 Nov), during which more than 10 5 images and spectra were obtained. 

Prior remote sensing showed that Hartley 2’s nucleus has an average radius 1/5 that 
of comet Tempel 1 ’s nucleus (J,6), yet it releases more gas per unit time at perihelion, 
even allowing for the smaller perihelion distance of Hartley 2 (1.059 vs. 1 .506 AU). This 
puts Hartley 2 in a different class of activity than Tempel 1 or any of the other comets 
visited by spacecraft (Fig. SI ). The two comets have very different surface topography 
(Fig. 1) but whether the different topography is related to the hyperactivity is still being 
investigated. 
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Figure 1 Comparison of a small part of Tempel 1 (left) with Hartley 2 (right) at approximately 
the same image scale and with nearly identical instruments. Left, Impactor Targeting Sensor 
(ITS) image iv9000675, 9.1 m/pixel. Right, MRI image mv5O04032, 8.5 m/pixel. Sun is to the 
right. 


The Nucleus 


Spin state and variations 

The rotation state of the nucleus distinguishes the morning from the evening 
tenninator and, by comparison with longer-term coma observations, allows the number 
and relative strengths of active areas to be determined. Knowledge of the nuclear spin can 
also put constraints on internal distribution of mass in the nucleus, internal energy 
dissipation, and the magnitude of the net torque. 

The large variations in brightness in Figure 2, reduced to a measure of the amount of 
dust leaving the nucleus, show a period of roughly 18 hours but the spacing of peaks in 
the light curve show's a clear pattern that repeats every three cycles. We interpret this 
(SOM) as an excited state of rotation with each cycle corresponding to precession of the 
long axis of the nucleus around the angular momentum vector, with a period of 1 8.34 
hours at encounter. The pattern of 3 cycles is due to an approximate commensurability 
between this precession and the roll around the long axis with a period of 27.79 hours 
(55.42 hours is also possible; the ambiguity does not affect any conclusions in this paper). 
The orientation of the angular momentum vector is not yet tightly constrained but is 
within 10° of being perpendicular to the long axis. This excited state also implies a 
nodding motion of the long axis relative to the angular momentum vector but the 
observed near-axial symmetry of the shape limits this to an amplitude <1 The 
precession period is increasing at 0,1% per period near perihelion, an unusually high but 
not unprecedented rate of change for a comet. The roll period is decreasing. These 
changes are presumed due to torques produced by the outgassing. 
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Figure 2. Variation of visible flux with time in a 191-km square aperture at the comet. Flux is 
proportional to the amount of dust leaving the nucleus if physical parameters of the dust do 
not change (no changes in dust color or other properties appear in our data nor are any 
reported by Earth-based observers). The dotted box is expanded later (Fig 7). Near closest 
approach, the flux from the nucleus is too large to allow accurate photometry of the coma. 

Nuclear Shape 

The nucleus is bi-lobed in shape (Fig. 1), with a maximum length of 2.33 km. The 
shape is well constrained by stereo viewing of nearly half the object and for much of the 
remainder sampled by silhouettes against the light scattered from, the coma (Table 1). 

This bi-lobed shape is crudely similar to that of comet Borrelly ( 18 ), but is both relatively 
and absolutely smoother. The rotation is slow enough that gravity is sufficient to hold the 
two lobes together for any bulk density >100 kg m \ 


Table 1, Properties o 

Hartley 2's Nucleus 

Volume 

0.82+0.08 km 

Mean Radius 

0.58±0.02 km 

Diameter 

0.69 to 233 km 

Cross Section 

0.43- 1.59 km 2 

Surface gravity 

0.0013-0.0033 cm s 2 

Uniform density 

220 kg m' 3 

Precession Period 

18.34±0.04 h (Nov 4) 

Roll Period 

27.79 or 55.42±0.1 hr 

Geometric Albedo 

4% (average) 

Albedo range 

4x 
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Constraints on Density 

The fast ( 12.3 km's) flyby did not permit determination of the nuclear mass from the 
spacecraft’s trajectory. The smooth shape of the “waist” region connecting the two lobes 
might indicate material collecting in a gravitational low, such as observed on asteroid 
(25143) Itokawa ( 19 ), Collection could proceed by in-falling material landing in the 
gravitational low and/or by in situ fluidization of regolith induced by outflowing gas ( 20 ). 
If some form of frictionless, fluidized flow is responsible for the formation or 
modification of this region it should represent a ‘flat’ surface such that it lies along an 
equipotential with respect to the combined forces of both gravity and rotation. Under 
these assumptions, which may not be valid, the density of the nucleus can be estimated 
by fitting potential contours to the observed geometry of the waist (SOM). 

We assumed internal homogeneity, a precession period of 1 S.34 hr, and a wide range 
of densities. The variance is minimized for a bulk density of 220 kg nf 3 (Fig. S4). Even 
this minimum residual leaves large-scale slopes of up to a few degrees relative to the 
equipotential. A reasonable lower limit for the density is 180 kg m 1 because the waist is 
no longer a gravitational low for lower densities. The upper limit is not well determined 
but a 4x increase above the best fit density to 880 kg m ’ requires modest porosity for 
pure ice and significant porosity for plausible rock/ice mixtures. 

Any interpretation rests on the character of the surface at the waist, which is mottled 
on horizontal scales of 10-30 m, and has some isolated cases of local relief > 10 m. The 
mottling, local topography, and generally gradational boundary distinguishes the waist 
from the best observed flow' on Tempel 1 (7 ) and the ponded materials on (433) Eros ( 26 ) 
and on Itokawa ( 19 ) at least locally, but does not rule out the possibility that the overall 
shape of the waist is approximately an equipotential. If the surface in this region 
approaches an equipotential and was formed by flow's or by deposition similar to those on 
other objects, it has subsequently evolved, suffering several m of erosional etching. 

These facts suggest that the equipotential assumption may not be reliable and that, unlike 
the case for Tempel 1 w'here a different and more direct approach was possible ( 21 ), the 
density might be significantly higher than deduced under our assumptions. 

Geology of Surface 

The nucleus has two primary terrain types (Fig. 3a): knobby terrain characterized by 
rounded to angular elevated forms up to 50 m high and 80 m wide and relatively smooth 
regions occupying both the waist (Fig. 3c) between the two lobes and parts of the larger 
lobe (Fig. 3b). Tire elevated forms appear to be the larger members of a population with 
most examples near or below our practical mapping resolution (-12 m; e.g. Fig. 3d). 

The smoother areas have darker central regions, are elongate, and are partially bounded 
by strings of the elevated forms that make up the knobby terrain. There is marginally 
resolved mottling of the smooth regions on the larger lobe, and better resolved mottling 
and local topography of 10-30 m scale in the waist. The darkest regions of the larger lobe 
are slightly more sharply bounded than is the darker band within the waist. The elevated 
forms that constitute the rough terrain arc in some areas aligned along boundaries of 
albedo markings, and follow' much of the southern edge of the “waist.” Many of the 
elevated forms exhibit 2 to 3 times higher albedo than the average, a much greater range 
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than seen on Tempel 1 . Ragged, somewhat sinuous, narrow depressions are visible at 
high incidence angles near the southern end, about 10 m deep, up to 90 m wide, and 
extending for over 250 m. 


Figure 3, Hartley 2, image mv6000002, ~ 7 m/pixel. Sun is to the right and the positive 
rotational pole is at the smaller (right) end. The view is from latitude ~ -3 3°. 


The average geometric albedo of the nucleus is ~- 4 %. Within the larger lobe area are 
several roughly equidimensional spots < SO m across that appear even darker than the 
larger, more elongated “dark'’ areas mentioned above (Fig. 3e). The darkest unshadowed 
spots are less than half the average brightness. While we can not rule out steep sided 
holes for the dark spots near the terminator, they generally occur in regions that are 
smooth in stereo and show no shadow signatures. Thus local albedos span at least a 
factor 4, compared to <2 on Tempel 1 . 


Jets occur in all terrains, but are clustered in the rough topography of the smaller lobe 
and mid- to northern part of the larger lobe (Fig. 3, also Fig. 5). Such clustering of jets 
has also been seen at comet 1 P/Halley ( 27 ). At least some jets appear to originate at or 
near large, bright elevated forms. Jets also originate beyond the terminator in areas with 
no direct sunlight, such as along the lower edge of the larger lobe in Fig. 4. Even our 
resolution of 1 0-12m is not sufficient to clearly resolve the morphology of the sources of 
the jets. 


This cornet lacks a population of depressions such as those that dominate 8 IP/ Wild 2 
( 28 ), or those scattered across Tempel 1 ( 1 . 29 ). The knobby terrain is similar to some of 
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the rougher areas on Tempel 1 . The smoother regions on Hartley 2 do not show the 
striations suggestive of flow markings on the best observed such region on Tempel 1 , and 
they are more gradationally bounded. The combination of terrains is very different from 
Tempel 1 or Wild 2, and this comet lacks exposures of thick, internal layers that were 
prominent on Tempel 1 . 


Nuclear Activity 

CN Anomaly 

Gaseous CN abundances were measured routinely from the start of observations on 5 
Sep (SOM). The long-term CN gas production gradually increased from 6xl0 24 s' 1 on 5 
Sep to a peak of 3 X lO" 5 s ! at perihelion (28 Oet), after which it decreased again to 
~2.4*I0 25 s *at closest approach and J.2*10 25 s' 1 on 25 Nov (Fig. S5). During most of 
the encounter the CN production varied periodically with the precession of the nucleus, 
as did the grains and other gases. 

The production of CN also exhibited an anomalous increase, by a factor of ~7, 
between 9 and 17 Sep. after which it slowly decreased, returning to the long-term trend 
line by 24 Sep (Fig. S5). There was a very weak increase in dust release (AA(0)fp <10 
cm) above its trend in that same period. There was no corresponding increase of water 
( 34 ). Integrating over the period 9-24 Sep. and subtracting a baseline gas production rate 
of7x!0- 4 s' 1 , we find that ~2xl0‘' 1 CN radicals (~800 tons) were released in the anomaly. 

The long-duration, gradual increase and decrease of gaseous emission without a 
corresponding increase in the dust production is atypical of cometary outbursts, which 
have sudden onsets and are usually accompanied by considerable dust. The increase is 
unlike the activity' observ ed at 9P/Tempel 1 or the behavior of any other comet. The 
spatial profile of CN during the anomaly was very different from that during the rest of 
the observations and suggests formation of CN from some extended source other than 
photodissociation of HCN. This could be grains too dark to scatter much sunlight such as 
HCN polymers ( 35 ) or the CHON grains found at 1 P/Halley ( 36 ). They would need to be 
lifted by something abundant and volatile. 

Large Chunks 

Radar observations in October showed an ensemble of particles > a few cm ( 12 , 37 ). 
While clouds of large particles had been reported previously from radar measurements of 
other comets (35), the location relative to the nucleus and the composition are not known. 
Previous searches with remote sensing for an icy grain halo in comets have rarely been 
successful with the only detections being at large heliocentric distances ( 39 - 41 ). 


At EPOXI’s close approach, individual chunks were seen near the nucleus in many 
images from both MR1 and HR1 (Fig. 4). The motion of the spacecraft allowed 
determination of the positions of individual chunks and their motions. A sample of 50 
chunks has been followed in many different MRl images and, other than one at 28 km. all 
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were found to be within 1 5 km of the large end of the nucleus. The motions are all very 
slow, 80% moving at less than 0.5 m s' ! and the fastest moving at <2 ms' 1 . This implies 
minimum lifetimes of 10 4 s for many of the chunks. Escape velocity is poorly defined 
near the surface because of the rotating, elongated shape as well as the uncertainty in the 
mass of the nucleus. We estimate that 10-20% of the chunks are moving at less than their 
local escape velocity. 

Sizes are estimated from the brightness of the chunks (SOM). The apparent flux from 
each of >10 4 individual chunks ranges from -10' B to <~1 O'” W m 2 pm' in the HR1 
image in Fig. 4, similar to that measured in the MR1 image taken close in time (Fig. 4). 
Sampling below 10 i2 W m " uni 1 is incomplete. We considered two extreme cases for 
the scattering properties - icy chunks scattering with the albedo and phase function of 
Europa (42-43) and dirty chunks scattering with the albedo and phase function of the 
nucleus of comet Tempel 1 (44). The range of measured fluxes corresponds to radii of 
10-150 cm if they are dirty chunks and 1 - 1 5 cm if they are ice (nearly pure). Below the 
minimum size the chunks blend into the background of unresolved, smaller chunks or 
grains. Since meter-sized objects are at the extreme of what can be lifted by gas drag and 
since the smaller, unresolved chunks are demonstrably icy (see below) we argue that the 
largest chunks are icy and roughly 10-20 cm in radius. 


The size distribution implied by the fluxes of the discrete chunks is unusually steep. 
Most of the mass and most of the cross-section are in the smallest grains. The discrete 
chunks contribute roughly 4% of the total surface brightness in the innermost coma (<5 
km) and those >~5 cm (completeness limit) are widely spaced at 4 X 10 * m \ If we 
extrapolate the size distribution (SOM) we find that >100% of the surface brightness is 
accounted for just with chunks > 0.5 mm. The total cross section of discrete chunks, 
roughly the same size chunks as ones to which radar is sensitive, is very' small compared 
to that detected by radar (12, 3 7), which presumably is detecting much darker chunks 
over a much larger field of view. 



Figure 4. Original HRI image (left, hv5004024, E-66s, range 915 km). Deconvolved image in 
center. Rightmost panel is an MR1 context image (mv.5004029) showing the location of the 
HRI field above the large lobe of the nucleus. Arrows indicate projected directions to Sun and 


Earth. 
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Heterogeneity of Dominant Volatiles 

Spectral scans of the comet were obtained from 1 Oct to 26 Nov, including several in 
which the nucleus was spatially resolved. Figure 5 is from a scan taken 7 minutes after 
closest approach with an MRI image taken at nearly the same time. Red boxes show 
regions where we have extracted the two spectra shown in Figure 6, both of which have 
had the continuum removed manually. The ratio of the FFO band to the CO? band varies 
by 2.9x spatially. In these maps (Fig. 5), the emission bands of HjO and CCF are both 
somewhat optically thick, implying that variations in column density could be larger than 
in brightness. 



Figure 5. Relative spatial distribution in the coma of Hartley 2. The red boxes (5x5 pixels; 52 
m/pixel) indicate regions sampled to produce the spectra in Figure 6. Panels labeled COz, 
Organics, and H 2 O Vapor are maps of the total flux in the relevant emission bands. The panel 
labeled HzO Ice is a map of the depth of the ice absorption feature at 3 pm. Each panel has 
been individually linearly stretched. All spectral images are from a scan at E+7 min, 
hi50060Q0, Sun is to the right. 
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Figure 6. Absolute spectra in the coma of Hartley 2 at CA+7 min. The black curve is the water- 
rich but ice-free spectrum from the box in the panel “H 2 O Vapor" in Figure 5, above the waist 
of the nucleus. The red curve is the ice-rich spectrum from the box in the panel "H 2 0 Ice” in 
Figure 5, the jet region near the small end of the nucleus. 

The maps show a water-vapor-rich region extending roughly perpendicular to the 
waist of the nucleus and presumably arising from the waist. This region has relatively 
little C(>, relatively little gaseous organics, and thus far no detectable water ice. The 
region of the jets off the end of the smaller lobe of the nucleus is rich in CO?, rich in 
organics, and rich in water ice but with a lower column density of water vapor than above 
the waist. There is significant ice in jets emanating from beyond the terminator along the 
lower edge of the larger lobe of the nucleus. The boundaries and the direct association 
with the major units of the nucteus seen here are dramatic and suggest very different 
histories for the waist and the remainder of the nucleus. The coincidence of strong 
absorption bands of ice with jets that are bright in the continuum suggests that jets are 
bright when highly reflective ice is present in the jets, and conversely that jets are usually 
fainter than the nucleus because they are of optically thin, relatively dark material. 

Theoretical calculations of scattering by icy grains (SOM) show that the predominant 
scattering grains must be smaller than 1 0pm. However. > 1 00% of the surface brightness 
can be accounted for by extrapolating the chunks to a size of 0.4 mm, implying that the 
chunks are fluffy aggregates or clusters of ~1 -pm solid grains. Either most of the 
aggregates of order 1 -mm have broken up or they mimic the scattering of the small 
grains. This result is very similar to the result obtained at Tempel 1 after the impact (no 
ice was observed prior to the impact). Those grains were predominantly micron-sized 
( 47 ). The similarity' between excavated material from Tempel 1 and ambient outgassing 
from Hartley 2 suggests that the constituent grains of solid iec are of order a micron in 
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most comets. Based on calculations of lifetimes ( 48 - 50 ) for the <10-pm solid 
components, the icc must be nearly pure for the grains to persist. 

The detection of strong absorption by ice, the detection of very large chunks in the 
coma, the concentration of all species other than H?0 vapor away from the waist of the 
nucleus, and the relatively smooth surface of the waist lead us to suggest that the material 
at the waist has been redeposited as a mixture of dirty grains and fluffy, icy aggregates 
that have not yet sublimed. The warmth of the dirty grains then leads to sublimation of 
the icy grains just below the surface. We conclude that this aspect of the chemical 
heterogeneity of the nucleus of Hartley 2 is likely evolutionary 7 . 

To determine the absolute abundance ratio, we considered a spectral map made three 
rotations (55 hours) earlier, when both the precession and roll orientations were the same. 
Spectra were extracted from a 1 20-km box and from a 600-km box, both centered on the 
brightest pixel of thermal emission (a better proxy for the nucleus than a reflected light 
center). In an aperture of 600*600 km centered on the nucleus (Fig. SB), and assuming an 
outflow speed of 0.5 km/s we find average production rates, QfyfyO) = 1 .0* 10 :s s’ 1 and 
Q(C02) = 2.0* 10 27 s' 1 for -20% fraction of CO 2 . This is higher than the fraction 
obtained in previous measurements of the global production of CO 2 in this comet ( 51 - 53 ). 

In Figure 7 we compare a portion of the visual light curve with the variation of CO 2 
and FTO from the spectral scans. The scale is arbitrary so only relative variations are 
meaningful. The ratio CCfyTfyO varies by a factor 2 between maxima and minima. The 
lower portion of Figure 7 shows images of the CO? and the H 2 O from the spectral maps. 
The red line indicates the position of the nucleus as defined by the peak thermal pixel. 
Close inspection shows that CO? is more sunward (up in the figure) than I 1?0 near the 
maxima, reflecting the different spatial distributions. This suggests that the CO 2 /H 2 O 
ratio is less in the large lobe of the nucleus than the small lobe, but this is a very 7 tentative 
conclusion until the rotational state is fully understood. If true, this heterogeneity is 
almost certainly primordial, unlike the ambiguous interpretation for the heterogeneity of 
Tempel 1 ( 52 ). 
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Figure 7. Light curve in visible light compared with light curve in H 2 0 and C0 2 , The curve for 
grains is an expanded version of the box in Figure 2. The data-points for gas are derived from 
spectral maps. The lower portion shows an expanded view with images of the C0 2 and H 2 0 for 
each point in the dotted box. The centroid of reflected continuum in each scan has been 
aligned with the red, horizontal line. The sun is up in the small images of gas. 


Summary & Conclusions 

Comet 103P/Hartley 2 differs in many ways from 9P/Tempel 1 and is an ideal 
example of hyperactive comets, ones that produce more H 2 0 per unit time than should be 
possible by sublimation from the small surface area of their nuclei. Super- volatiles, 
specifically COi in the case of Hartley 2, are the primary drivers of activity’. The super- 
volatiles drag out chunks of nearly pure water-ice, which then sublime to provide a large 
fraction of the total i frO gaseous output of the comet. Other hyperacti ve comets include 
46P/Wirtanen and 21P/Giacobini-Zinner. 

In Hartley 2. H 2 O sublimes from the waist with much lower content of C0 2 and 
barely any trace of icy grains. We tentatively interpret the waist as a secondary' deposit 
of material, although the mechanism of redeposition remains unclear. The most likely 
mechanism involves fallback of both refractories and icy chunks from the periphery of 
the active regions. C0 2 /H 2 0 varies by a factor 2, likely between one end and the other. 

From HST measurements (55), CO < 0.3%, implying CO 2 CO 60. a far more 
oxidized environment than is consistent with any model of the outer protoplanetary disk. 
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This is also different from Tempel 1 (~1) (54) and Halley (<1 ) (56). A large anomaly in 
CN, too slow to be called an outburst, is unexplained. 
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Supplementary On-line Material 


Global Comparison 



Figure SI. Images at the same scale for all cometary nuclei observed by spacecraft. 
Differences in overall shape are dramatic, as are the differences in contrast between the 
nuclei and their associated jets, which are brighter than the nucleus at Halley and Hartley 2 
and much fainter than the nucleus at the other comets. 
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Rotational State 


For background discussion of excited states of rotation in the context of small bodies, 
and in particular the choice of terminology, see (7-9), These references also include the 
background mathematical derivations. Our discussion is based on a Long Axis Mode 
(LAM) as described in those references since this is far more likely dynamically possible 
than a short-axis mode (SAM) based on the observed properties, including the observed 
axial ratios. 

Periodogram analysis of the light curve on approach revealed seven periods in the 
earliest data, from E-60 to E-40 d (E ~ time of encounter), of which six are shown in 
Figure S2. Subsequent accumulating data then showed that these periods were changing 
over the course of the encounter (Fig. S2). Periods P2 and P4 increase with time and are 
associated with the precession period of the long axis around the angular momentum 
vector and its first sub-harmonic. Thus the period of precession of the long-axis is 
varying from <17 to >18 hours, with a value 1 8.34±0.04 h at encounter. The observed 
motion of the nucleus in the short interval over which the nucleus was clearly spatially 
resolved is consistent with P2 but not with P4 as the true precession period. Furthermore, 
the morphology of the coma as seen from ground-based telescopes also repeats 
approximately (but not exactly) with the period P2. The variation of P2 with time is 
shown in more detail in Figure S3, w'here we have included several published values 
from Earth-based observatories. 

Periods PI, P3, P5, and P6 decrease with time (Fig. S2, Table SI) and are associated 
with tire roll around the long axis. We adopt P3 = 27.79±0.31 h at encounter as the roll 
period. PI has very low power in the periodogram analysis and is therefore likely a 
harmonic, but P5 has power comparable to P3 and can not be excluded as the true period. 
Since the roll around the long axis has low linear speeds at the surface, the motion can 
not be seen in our resolved images. The ambiguity has no effect on any of our 
conclusions in this paper. The rates of change of the periods, determined from fitting a 
straight line to the data in Figure S2, are given in Table SI . P? is relatively weak, 
approximately 8*P2 and off scale in Fig. S2. 

As a first approximation, the spatial orientation of the rotation state has been 
determined assuming that the end of the long axis traces out a great circle on the sky. 

This is equivalent to assuming that the rotation is not excited, a reasonable first 
approximation since our subsequent studies show that the excitation is small. We 
constrained the direction of the total angular momentum vector by the change in position 
angle of the projected axis of figure of the nucleus and its apparent length on approach to 
encounter. In a set of 1 1 resolved and deconvolved HR1 images 1 from the last hour of 
approach, the actual orientation of the long axis moved over an 1 8° arc on the sky, which 
constrained the direction of the total Angular Momentum Vector (AMV) to RA, Dec = 
59°, 55°± 2° in this approximation. At encounter (JD2455505.G8318) the small, active 
end of the long axis points towards RA, Dec = 226. 1 39.4° (± 0.7° circular error). 

Although the roll period is still uncertain it is possible to conclude that the excitation 
is not extreme based on ratios of the principal moments of inertia obtained from the shape 
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model under the assumption of a homogeneous mass distribution in the interior. Using 
the theory of force-free rigid body rotation ( 8 ), we find that the long axis should be 
inclined to the AMV by >80° and that any “nodding” motion of this axis should have an 
amplitude < 1°, too small to be detectable in our data. The excited spin state is still being 
evaluated and we expect that the excited solution will displace the orientation of the 
AMV by <~20°. 



Figure S2. Variation of the periodicities seen in the light-curve of the comet through the 
EPOXI/BIXI encounter. The periodicities P2 and P4 increase with time {a third, longer period 
that increases with time, P7, is not shown here in order to increase the clarity of the 
presentation); the remaining periodicities, PI, P3, PS, and P6, decrease with time. We 
associate the former group with the precession of the long axis of the nucleus around the total 
angular momentum vector and the latter group with the spin component around the long axis. 
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Weighted Trend-Line Parameters for 103P/Hartley 2 


Frequency 

Period 

# 

Si ope x 10 6 (hr ydy) 

Intercept (hr*) 

Slope x 1(P (hr ydv} 

Intercept (hr) 

fi 

+2 ±17 

0.0901 + 0.0004 

-0.5 ± 2 

11.09 ±0.05 

h 

-72 + 7 

0.0543 ± 0.0001 

+24 + 2 

18.34 ± 0.04 

k 

+ 21 + 17 

0.0360 + 0.0004 

-17± 13 

27.79 + 0.31 

f4 

-37 z 3 

0.0272 + 0.0001 

+49 i 4 

36.73 ± 0.08 

fs. 

+4 + 2 

0.0180 1 0.00003 

-5 ±7 

55.42 + 0.10 

f<, 

+22 ±S 

0.0121 ±0.0001 

-165 ±32 

82.64 ± 0.57 


Table SI. Weighted slopes and S2 and S3. The intercepts represent our best estimates of the 
component periodicities at encounter. The form of the trend-line is y = slope • x+ yo where x is 
days from encounter and y 0 is the intercept at encounter (JD 2455505.083 18).The number (#) 
of the periodicity is given in column 1 and shown in Figures S2 and S3. 

Our estimates of the total angular momentum and total rotational energy per unit 
mass are -1 . 1 x 10 5 m 2 h 1 and ~2>' 10 4 nrh' and these are changing at the rather 
substantial rates of -0. 1 %/d and ~0.2 %/d. respectively. The change in rotational period 
(or angular momentum) of Tempel 1, for comparison, is -0.02% per day near perihelion 
( 10 ). The harmonic relationships between the periodicities are quite clear in the data with 
a strong first sub-harmonic of both the precession and roll periods appearing in the 
periodograms. The rapid change in the periods and the temporal changes in the non- 
gravifational acceleration of the comet in its orbit both are indicative of the combination 
of a high outgassing rate coupled with a small mass and, for the periods, a very elongated 
shape. Figure S3 combines our data on the preeessional period, P2, with a variety of 
other observations ( 11 - 15 ) and shows that the preeessional period remains constant while 
the comet is far from the sun and inactive but that it changes rapidly around perihelion 
while the comet is active. Finally, the precession and roll periods appear to be nearly 
commensurate (-2:3), but this is likely a coincidence since these periods are changing 
rapidly in opposite directions. Such commensurability has been seen before in the case of 
1 P/Halley ( 8 , 16 ), where it is apparently a stable resonant phenomenon. It is a feature 
that has been seen in numerical simulations of the evolution of cometary spin ( 1 7). 
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Figure S3. Change in the precessional period of the long axis around the total angular 
momentum vector as seen from Earth and from the Deep Impact spacecraft. Color code for 
Earth-based observations: grey (17); aqua (12); red (13); gold (14); green (15). The seven 
EPOXI/DIXI points are in black. 


1 HRI Images used were hvS000005, hv5000022, hv5000065, hv5000095, 
hv5000118, lw5000139, hv5000152, hv5000l52, hv5002000, hv5002015, 
hv5002028, and hv5002044. 
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Equipotential Surface 

Thc equipotential model for the surface gravity was used to fit the shape of only the 
waist. We have not attempted to model the shape of the entire nucleus as an 
equipotential. The model for the potential includes the effect of the precession of the long 
axis around the Angular Momentum Vector but does not include the effect of the roll 
around the long axis. Because the period of roll is long and the lever-arm at the surface is 
short, the contribution of roll to the potential is unimportant. Figure S4 shows the 
normalized variance' between the observed shape and the calculated potential, in the 
waist only, as a function of density. The minimum variance is at p = 220 kg uf 3 , which 
corresponds to a total mass = 1.84x 10 11 kg. At that minimum, the normalized variance 
indicates 1 1% variations of the actual potential around the equipotential. As shown in the 
figure, the fit deteriorates rapidly for bulk densities less than the best fit value and the 
waist ceases to be a gravitational low for bulk density <180 kg m which would 
invalidate the assumptions. This sets a lower limit on the density if our assumption of an 
equipotential is correct. The two lobes wrnuld be gravitationally bound to each other, with 
no strength required, even for densities down to 1 00 kg m‘\ The fit degrades much more 
slowly for higher values of density. Given the deviations from an equipotential, densities 
>220 kg nr 3 are plausible but even 4x higher densities require moderate to high porosity 
of the bulk material for reasonable assumptions about the rock/ice ratio of the nucleus. 
The density is clearly in the range found for other, well studied comets such as Tempel 1 
at 400 or 450 kg m' 3 ( 21 , 22 ), 81P/Wild 2 at 600 kg m" 3 ( 23 ), and 19P/Borrelly at 240-490 
kg m 3 ( 24 , 25). All these determinations are model dependent and have large 
uncertainties but they are all consistent in being much less than the bulk density of pure 
ice. 
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Figure S4. Variance of the potential fitted to the waist of Hartley 2 as a function of assumed 
density. 
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CN Anomaly 

Routine imaging of the comet began on 5 Sep ~ E-60d, using primarily a clear filter 
sensitive to reflected light from grains in the coma and a narrow-band filter isolating an 
emission band of gaseous CN at 387 nm. The narrow band images were used to derive 
absolute gas and dust production rates of Hartley 2. To remove the continuum 
contribution in the CN filter, we assumed an unreddened solar continuum and convolved 
this with the filter transmission and quantum efficiency of the detector ( 4 ). The resulting 
CN fluxes were used to derive column densities using fluorescent efficiencies scaled to 
the comet’s heliocentric distance and velocity ( 30 ). To derive gas production rates, we 
compared the measured column densities to an image created from a Haser model ( 31 - 
32 ) with parent scale length 2*]0 4 km & daughter scale length 2x10" km, values 
consistent with the expected lifetimes of HCN and CN against photodissociation. While 
the spatial distribution of CN released during the comet’s regular activity can be 
described well with this Haser model, the spatial distribution of CN during the anomaly 
has a significantly shallower slope, suggesting an additional, extended source. A 
comparison of production of CN (Fig. S5) with the production of water ( 34 ) indicates that 
there was no corresponding increase in water release. A large amount of CN with little or 
no corresponding dust or water is quite surprising. The spatial profile of CN during the 
anomaly suggests that the CN was in grains too small or too dark to scatter much sunlight 
such as HCN polymers (35) or the CHON grains found at 1 P/Halley ( 36 ). They would 
need to be lifted by something abundant and volatile. We do not yet have an appropriate 
physical model for the anomaly since the mechanism for producing the CN in the 
anomaly is apparently different from the one that operates normally. 
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Figure S5. Temporal variation of CN and grain production. CN production in 
molecules/s compared with the production of grains measured bv the parameter Afp 
( 33 ). 
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Discrete Chunks 


The detailed size distributions and the spatial distributions of the large chunks are 
shown in figure S6 for a particular MR1 image, mv5006001, taken ~1.5 minutes after 
closest approach. Fluxes were measured in several images in three ways: 1) by aperture 
photometry of each chunk in calibrated HR1 images, 2) by aperture photometry of each 
chunk in deconvolved, calibrated HRI images, and 3) by fitting a point-spread function to 
individual chunks in calibrated MRI images. All three approaches lead to consistent 
results. Completeness was assessed by randomly inserting artificial sources of varying 
brightness and then finding them in the images using our automated routines. In the size 
distributions in Figure S6 (upper two panels), points are plotted both as directly measured 
and after being corrected for incompleteness in recognizing chunks with our software. 
The same counting data are fitted separately to number as a function of measured flux 
(proportional to the cross-section of the chunk) and to number as a function of diameter 
(assuming Europa-like, icy scattering functions: 42 - 43 ). We also tried dirty scattering 
functions based on the nucleus of comet Tempel 1 ( 44 ) but consider that less likely as 
discussed in the main text. Note that the chunks are distributed all around the nucleus 
without any prominent signs of discrete sources on the nucleus, suggesting that the 
motions of the grains have become isotropic in the nuclear frame due to the frequently 
varying gravitational field with every precession period. 
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Figure S6. Distribution of chunks seen in MRI image mv5006001. The image was taken 1.5 
min after closest approach. Upper left shows the distribution as a function of measured flux 
from each chunk ]W nr 2 pm 1 ] while the upper right shows the distribution as a function of 
diameter [cm] assuming icy chunks. The lower left is an empirical test of the completeness of 
our counting as a function of the observed flux and is used to correct the actual data in the 
upper plots. The lower right shows a map of the spatial distribution of the particles in this 
image with sun at the right. 

With an expansion velocity upper-limit of 30 cm s' ! , we compute a mass-loss rate in 
the large chunks alone of 0.5 to 2 x 10" kg s~ ! . but the uncertainties due to possible 
sublimation of icy grains and the unknown density of aggregates could alter this number 
substantially. 
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Spectral Maps 

Spectral processing for each long-slit spectrum includes: a) removal of instrumental 
thermal background that depends strongly on known quantities such as temperatures 
inside the spectrometer but also when the most recent previous spectrum was taken and 
sometimes on the signal in the last previous spectrum; b) correction for non-linearity in 
the gain of the spectrometer; c) conversion to radiance. Extraction of spectra includes 
allowance for curvature of the slit as imaged onto the detector (commonly known as 
spectral “smile”), summing or averaging over spatial pixels both along the slit and from 
one spectrum to another during a scan. For most purposes, more useful spectra are then 
obtained by subtracting the continuum, which consists of reflected sunlight at the shorter 
wavelengths and thermal emission at the longer wavelengths. For the spectra presented 
in this paper, the continuum has been removed by linearly interpolating between 
segments of the spectrum thought to be representative of the continuum (reflected 
sunlight). When only absorption features are of interest, the spectrum is divided by the 
solar spectrum before subtracting the continuum but this step is omitted when the fluxes 
of emission features are also of interest. Eventually a proper thermal and scattering 
model will be used to remove the continuum but this is a major task. From past 
experience we have limited our conclusions here to ones that will not be affected by 
improvements to the processing. 

Figure 5 was from a spectral scan centered at E+7min consisting of 56 spectra taken 
over ~2.5 minutes with an exposure time of 1 .4 sec per frame during a continuous scan at 
1 slit-width per frame. The scan will be available in the PDS archive as image id 
165006000. The maps shown are integrated over the relevant spectral range. A mask was 
used to suppress the nucleus. The maps have a scale of 52 m/pixel; the maps of the gases 
are smoothed with a 3x3-pixel box but the map of the ice is not smoothed. All maps are 
linear in intensity but with varying amounts of stretch. The stretch in the images is best 
understood by looking at the actual spectra in Figure 6. Those spectra are averages of the 
surface brightness (radiance) over the boxes and were reduced from the calibrated 
radiance spectra by fitting straight lines to portions of the spectrum thought to be true 
continuum. Radiances integrated over the emission bands in the two boxes are: Water- 
rich waist: B(H 2 0)= l,7*10'-\ B(C0 2 ) = 4.1 xlO" 4 , and Ice-rich lobe: B(H 2 0) = 1 . 1 * 1 0 \ 
B(C0 2 ) = 8.1 x I CP, all in [W nf" sr' 1 ]. At this stage of processing, the spectral profile of 
the organic emission is not yet reliable due to difficulty in separating the ice absorption, 
the smooth reflected continuum, and the thermal emission, all of which matter at the 
wavelength of the organics. Because these factors vary from pixel to pixel, the 
uncertainties in the map of organics are much larger than the uncertainties in the maps of 
the strong emission features of 1TO and CCX Based on abundances reported from 
remote sensing for this and other comets, we presume that the emission is dominated by 
gaseous methanol (CHjOH), but further discussion of the organic emission is premature. 

In order to study the icy grains further, we have calculated models for the reflection 
spectrum using the optical constants of ice (45) and compared them with the relative 
reflectance spectrum in Figure S7. The observed spectrum is the same as the one for the 
ice-rich region in Figure 6, but processed to relative reflectance. In this case, the thermal 
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emission was subtracted from the observed spectrum as an empirically fit Planck function 
and the result was divided by the solar spectrum so that any residual “color” due to the 
scattering properties would remain. It was then normalized at 1 .8 pm. The model 
spectra {46) are for solid spheres of pure, hexagonal ice in an areal mixture with particles 
assumed to have a reddish slope, characteristic of refractory grains in comets, i. e . , the 
sum of a reddish continuum and a theoretical ice spectrum. The best fit is sensitive to the 
sizes of the particles but a range of mixing ratios and slopes for the reddish material 
provides nearly equivalent fits. Archival ID for the data is hi5006000. 




Figure S7 Relative reflectance spectrum of the coma. The smooth curves are theoretical 
spectra of icy grains of 1 (red), 1 0 (green), and 100 (purple) pm diameter, normalized to unity 
at 1.8 pm. The latter two curves have been displaced vertically for clarity. The variation in 
relative strength of the absorption features at l.S, 2.0, and 3.2 pm shows that the dominant 
size of the ice must be less than a few pm. 

'The fits to the relative reflectance spectrum are primarily sensitive to the size of the 
icy particles and it is clear that the dominant source of scattered light from the ice must 
be particles with diameters less than 10pm. This is consistent with the earlier conclusion 
that the dominant source of light must be the smallest particles, no matter what the large 
chunks are. We note that large fluffy aggregates can mimic individual grains of the same 
size as the constituent pieces of the aggregate and can thus appear small in this type of 
analysis, even though the actual, porous particles may be much larger. This result is very' 
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similar to the result obtained at Tempcl 1 after the impact (no ice was observed prior to 
the impact). Those grains, which had undergone an excavating shock, were 
predominantly micron-sized (47). The similarity between excavated material from 
Tempei 1 and ambient outgassing from Hartley 2 suggests that the constituent solid 
grains of ice are of order a micron in most comets and that the icy grains seen at Hartley 
2 are also from the subsurface. 

The spectral map made at E-55 hr was taken from spectral scan hi4000039, integrated 
over areas of 120x120 km and 600x600 km, centered on the nuclear position (defined by 
the peak thermal emission) to produce the spectra in Figure S8. The continuum was 
determined manually and interpolated with straight lines. The interpolation uncertainties 
could allow some absorption by ice in the 3 -pm band, which would imply a higher total 
production of water vapor than deduced from these spectra. The band ratio of CCT/ITO 
differs between the two spectra indicating that the optical depth of the bands is still 
affecting the result in the 120-km box. Average column densities based on an optically 
thin assumption are: 120-km: N(H 2 0) = 6.2*10 t5 cm' 2 , N(C0 2 )= 8.4xl0 14 cm' 2 ; 600- 
km: N(H 2 0) = 1.7*10 15 cm' 2 , N(C0 2 ) ~ 3.2* 10 14 cm' 2 . For simple, radial outflow', the 
average column density should be inversely proportional to the size of box, but this can 
be altered either by an extended source of H 2 0 from the icy grains or by significant 
optical depth of either species in the 120-km box. For assessing the total release of 
volatiles, the 600-km box is the better choice. 



Wavelength (pm) 


Figure S8. Continuum-removed spectra of Hartley 2 taken 55 hours before closest approach. 

The spectra used in Figure 10 were from a large number of short, spatial scans taken 
at a rapid cadence on departure, hi4300007 through hi4500017. To ensure coverage of 
the comet by short scans, the rate was set to two slit widths per frame rather than the 
more usual one slit width per frame. Thus a square box of 5*5 pixels was used to extract 
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the flux but this actually corresponds to a 10x5 pixel field of view. The range was 
changing during the period of the scans so the size of the box is not constant in km but is 
roughly 350> : 175 km. large enough that the effects of optical depth are minimized. 
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